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Perturbation Theory for Fat-link Fermion Actions 
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We discuss weak coupling perturbation theory for lattice actions in which the fermions couple to smeared 
gauge links. The normally large integrals that appear in lattice perturbation theory are drastically reduced. Even 
without detailed calculation, it is easy to determine to good accuracy the scale of the logarithms that appear 
in cases where an anomalous dimension is present. We describe several 1-loop examples for fat-link Wilson and 
clover fermions. including the additive mass shift, the relation between the lattice and MS quark masses, and 
the ajcial current renormalization factor {Za) for light-light and static-light currents. 
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Smearing or "fattening" the gauge links in a 
fermion action has several advantages. Indeed, 
with sufficient fattening: 

• Exceptional configurations are suppressed. 

• Additive mass renormalization is small. 

• Finite matching constants (Z factors) are 
very close to 1. 

• We expect that tree level 0[a) improvement 
(clover coefficient Csw = 1) is close to all- 
orders (in g) 0{a) improvement. 

The last three features are closely related and 
can be understood in perturbation theory. Fur- 
thermore, these actions are already being used in 
Monte-Carlo computations {e.g., by MILC (|]), so 
it is important, as a first step, to have perturba- 
tive evaluations of the renormalization constants. 
Non-perturbative evaluations, where possible, are 
of course a very important next step. 

In the present work, we consider the ordinary 
Wilson and clover fermion actions in which the 
gauge links are fattened by APE smearing 
The link after m + 1 smearings is related to the 
link after m smearings by 



(m+l) 



{x)=v({l-c)vt\x) + 



+ Uj^™^ {x - z>)Fi™^ {x - i>)uj") {x-i) + fi)]) 

where V is the projection back into SU(3) and 
V^(x) = Ufi{x), the original ("thin") link. 

If one fixes c and performs a fixed number of 
smearings, A^, as a — > 0, the action remains lo- 
cal and the continuum limit will be unaffected. 
In practice, MILC has been trying c = 0.45 and 
A^ = 10 for heavy-light decay constants and form 
factors; we use these values for all "fat" results 
quoted in this paper. (Unless specified otherwise, 
we take Csw = 1 in clover computations.) 

We are interested in computing 2- and 4-quark 
operator renormalization/matching constants at 
one loop. For such computations, only the linear 
part of the relation between fat and thin links is 



relevant {Vi'i"^\x) 



exp(iav4^™^(a;)): 
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Quadratic terms in (^, which would be relevant 
for tadpole graphs only, appear as commutators 
and therefore do not contribute, since tadpoles 
are symmetric in the two gluons. 

In momentum space, eq. (hJ) becomes 
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where 

V(9) = /('z)(5..-%^) + %|^ 



(2) 



(3) 



2 



with = f sin(^) and /(g) = 1 - ^f. ^ 

After N smearings, h^v{q) becomes h'^^{q), 
which is just /i^^ with / replaced by f^. If 
< c < 0.75, then |/| < 1 over the entire Bril- 
louin zone, and the factor will, for large N , 
strongly suppress values of q outside of a small 
ball around q — Q. 

For small c, it is easy to find the effective range 
of the fattening: 

r(g)^exp(-^a^(-^)-^a^ (4) 

Thus, even an = 10, c = 0.45 fat link has a 
mean size of only a couple of lattice spacings. 

The longitudinal part of A^^ in eq. @j is unaf- 
fected by smearing, because it is a pure gauge, for 
which parallel transport is independent of path. 

The Feynman rules are easily derived from 
eqs. (|]Jq). Each quark-gluon vertex gets a form 
factor h^iyiq), where q is the gluon momentum. 
If all gluon lines start and end on fermion lines, 
then, effectively, the gluon propagator changes by 
Dfj^iy — > h^^Dxah^iy- The form of shows that 
Landau gauge is the most natural gauge for fat- 
link perturbation theory, since the Landau prop- 
agator will kill the longitudinal part of h. 

In many cases it is now easy to convert ordinary 
thin-link perturbation theory to the fat-link case. 
In particular, consider some pure lattice quantity 
{e.g., the additive mass shift, or Za for two light 
quarks) computed for thin links at one loop. The 
answer is of the form J ^^^^4 Iiattil)- If the thin- 
link calculation was done in Landau gauge with 
the loop momentum along the gluon line, then 
the fat-link result is obtained simply by replacing 

However, if the thin-link calculation was done 
in Feynman gauge with the loop momentum along 
the gluon line, the fat-link Landau gauge compu- 
tation will have some new terms, coming from the 
fN%^ part of hf^^iq). When / = 1, the sum of 
such terms must vanish by gauge invariance. If 
the cancellation occurs at the level of the inte- 
grands, multiplying by /■^^ will not affect it, and 
the change from thin to fat links is again accom- 
plished by Iiattiq) Iiatt{q)P'^{q)- In the cases 
we have checked (additive mass shift, multiplica- 
tive mass renormalization) , the cancellation does 



indeed occur at the level of the integrands. How- 
ever, we know of no general proof that this must 
be the case. If the cancellation required integra- 
tion by parts, for example, the fat-link computa- 
tion would contain some new terms not present 
in the thin-link case. 

We can now perform the lattice integrals nu- 
merically. An important first example is the sim- 
ple gluon tadpole {Iiatt = We find that it has 
the value 0.1547 in the thin-link case, but only 
0.0044 for fat links. The form factor /^^ has 
completely suppressed the UV part of the Bril- 
louin zone, which gives the dominant contribution 
to the tadpole. Thus tadpole improvement is un- 
necessary for fat links, at least with this much 
smearing. 

Another example is the axial current 

renormalization constant for two massless quarks. 
Writing Z'\ = l + ^Ca, we have 

r -15.80, thin Wilson 
(]i = J -13.80, thin clover (5) 
1-0.24, fat clover 

(Here clover currents are local: improved by the 
factor 1 -I- uiQ = 1 + ^ — 2^ ' ^'^^ derivative 
terms.) Again, the form factor has suppressed 
most of the integration region, making the inte- 
gral small. With this much fattening, all such 
pure lattice integrals should be highly suppressed. 
This explains why the additive mass renormal- 
ization is so small, and why a clover coefficient 
Csw = 1 should be very close to the nonpertur- 
batively improved value — certainly the tadpole- 
improved Csw is extremely close to 1. 

The matching of a divergent operator {i.e., one 
with an anomalous dimension) is slightly more 
complicated. In this case, a standard way to ex- 
press the lattice result in the thin case is: 

It^illatt {q)~ leant (?, 0)] 

+ J^Icontiq,aX), (6) 

where A stands for some IR cutoff, such as a quark 
mass or momentum, or a gluon mass inserted by 
hand. Icont{q,aX) is a simple continuum- like in- 
tegrand which has the same IR behavior as the 
lattice integrand, but which can be integrated an- 
alytically, giving the explicit In(aA) behavior. 
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In the fat case, the analogous formula is 

/(Cf^/'^ [hatt {q)- leant [qM 
+ I ■(2^Icont{q, a A) 
+ Jj^Icontiq,0){P^-l) , (7) 

where we have added and subtracted the integral 
of f'^^ leant- This separates the result into three 
terms: (1) a complicated, regularized lattice in- 
tegral, which is however suppressed for sufficient 
fattening and is thus numerically small, (2) the 
same continuum-like integral that appeared in the 
thin case, and (3) a simple integral involving the 
form factor, whose coefficient is determined by 
the anomalous dimension of the operator. Term 
(3) is not particularly small, since /^^ does not 
multiply the entire integrand. 

An example is the matching of the lattice mass 
to the continuum MS mass: 

2 

ZM(aM) = 1 + ^ (Cm - 6 log (af,)) . (8) 
We get 

[-12.95, thin Wilson 
Cm = < 19.31, thin clover (9) 
I -4.92, fat clover . 

Z^, the axial current renormalization constant 
in the static-light case, also has an anomalous 
dimension. It does not separate as easily as 
eq. (^. However, the calculations, following [Q, 
are straightforward. We find 

ZiiaMB) = 1 + ^ (C^' + 3 log (aMs)) (10) 

r -22.36, thin Wilson 
(^1 ^ } -16.41, thin clover (11) 
I 0.393, fat clover . 

(The clover currents are local.) 

The bad news about fat links is that, at least 
with the amount of fattening used in the examples 
above (TV = 10, c = 0.45), the lattice integrands 
are so suppressed at large q that they may become 
IR sensitive. A measure of the effective cutoff of 
the integrals is given by g* ||]. With Csw = Ij 
we find q*a = 0.71 for Z^; while q*a = 1.05 for 
So (which gives the shift in Kc). 



More direct evidence that perturbation theory 
for Eq is breaking down comes from a comparison 
with numerical results from the MILC collabora- 
tion. For /? = 5.6, m = .01, Csw = 1-0, N = 10, 
c = 0.45, the numerical result is Kc = 0.1256(1); 
whereas 1-loop perturbation theory, with q*a = 
1.05, gives Kc = 0.1251. Thus perturbation the- 
ory fails by a factor of six for Kc — 1/8. There are 
however sensitive cancellations at 1-loop making 
So anomalously small. 

For the same data set, Z]^ = 0.99 in perturba- 
tion theory. We guess this will be better behaved 
than So, because there are no delicate cancella- 
tions here. But even if — 1 is wrong by another 
factor of 6, the associated error in the final answer 
is only 5%, since Z^' is so close to 1. 

However, in cases with anomalous dimensions 
{e.g., Jb), the perturbative corrections are not 
necessarily small, and a failure of perturbation 
theory could lead to large errors. Some non- 
perturbative calculations arc needed to clarify 
this situation. 

In retrospect, somewhat less fattening in the 
MILC running would probably have been prefer- 
able. For example, for the clover Z^, N = 7, 
c = 0.25 gives q*a ~ 1.34 instead of 0.71, and may 
still suppress exceptional configurations enough. 
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